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ABSTRACT 

We investigate the relationship between brightest cluster galaxies (BCGs) and their host clusters 
using a sample of nearby galaxy clusters from the Representative XMM— Newton Cluster Structure 
Survey (REXCESS). The sample was imaged with the Southern Observatory for Astrophysical Re- 
search (SOAR) in R band to investigate the mass of the old stellar population. Using a metric 
radius of 12 h" 1 kpc, we found that the BCG luminosity depends weakly on overall cluster mass as 
Lbcg oc M° ; ' 18 °' 07 , consistent with previous work. We found that 90% of the BCGs are located 
within 0.035 r 50 o of the peak of the X-ray emission, including all of the cool core (CC) clusters. We 
also found an unexpected correlation between the BCG metric luminosity and the core gas density 
for non-cool core (non-CC) clusters, following a power law of n e oc L^cg ' 4 (where n e is measured 
at 0.008 ?"5oo). The correlation is not easily explained by star formation (which is weak in non-CC 
clusters) or overall cluster mass (which is not correlated with core gas density). The trend persists 
even when the BCG is not located near the peak of the X-ray emission, so proximity is not necessary. 
We suggest that, for non-CC clusters, this correlation implies that the same process that sets the 
central entropy of the cluster gas also determines the central stellar density of the BCG, and that this 
underlying physical process is likely to be mergers. 
Subject headings: X-rays:galaxies:clusters - cooling flows - galaxies:cD 



1. INTRODUCTION 

The intracluster gas in some galaxy clusters shows a 
central concentration in a cool core. The higher density 
of this core gas allows more rapid energy loss in the form 
of X-ray emission. In the absence of other energy sources, 
the classic "cooling flow" model suggests that the central 
gas should cool. The cool gas would rapidly condense and 
form stars at rates greater than lOOMQvr , but opti- 
cal colors and spectral lines indicate much lower rates 
of les s than ^10M Q yr _1 (reviewed in iDonahue fc Voitl 
12004 ). Thus, the gas must be heated by other processes, 
such as AGN activity or other forms of feedback that 
regulate the thermal prope rties of the gas (reviewed in 
iMcNamara fc Nulsenll2007ft . Yet the identity of the feed- 
back process is still widely debated, and current simula- 
tions tend to overpredict th e fraction of clusters with cool 
cores (e.o. iKav et al.l 12007( 1 . 

The gas core is often located near the Bright- 
est Clust er Galaxy (BCG ) , par ti cularly in cool core 
clusters (iJones fc Formanl 11983: IBildfell et all |2008; 
IRaffertv et alJ 120081 : ISanderson et all 12009ft . The BCG 
is not only luminous and massive, but has a more ex- 
tended optical light profile than other ellipticals, due to 
its unique merger history a t the bottom of the gravi- 
tational potential well {e.g. lHausman fc Ostrikcri ll978l : 
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I Vale fc Ostrikerl 12008( 1. Models predict that the BCG 
mass is correlated with the total cluster mass {e.g. 
iDe Lucia~fc Blaizot 20 Q7]), but the observed correlation 
is typically weaker (|Lin fc Mohrl l2004t iPopesso et al. | 
2007t |Brough et al.ll2008l : iWhilev et al.ll2008l ; lYang et al.1 
2008; iMittal et al.ll2009ft . 



To better understand the connection between BCGs 
and cool cores, we investigated the BCGs in 31 
southern hemisphere clusters from the Representative 
XMM- Newton Cluste r Structure Survey (REXCESS) 
(|Bohringer et al.ll2007f) . The REXCESS sample was cho- 
sen to evenly represent the range of X-ray luminosi- 
ties in the local (0.06 < z < 0.18) population, with 
no bias regarding X-ray morphology or central surface 
brightness. Thus, it provides an ideal laboratory for 
testing t he connections between BCGs an d their host 
clusters. IPratt et alJ (|2007H . iCroston et all (J2008f l. and 
iPratt et al.l ~(|2009) measure the X-ray luminosity, tem- 
perature, mass, core gas density, and cooling time. Here 
we present ground-based CCD (optical R band) imag- 
ing of the cluster BCGs and calculate colors relative to 
2MASS K band magnitudes. We investigate the connec- 
tions between the old stellar population and the proper- 
ties of the X-ray gas. 

Interestingly, the strongest correlation we found was 
not between the BCG mass and the total cluster mass, 
but between the central BCG stellar density and the core 
gas density. We describe the cluster sample and X-ray 
measurements in and the R band observations and 
BCG magnitude calculations in fj3] In ij4] and [5] we ex- 
amine this new correlation, as well as other connections 
between BCGs and their host clusters. Our conclusions 
are summarized in We use h — 0.7, fi m = 0.3, and 
£1a = 0.7 throughout. 

2. X-RAY SAMPLE 
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The REXCESS sample ()Bohringer et all 12007ft con- 
tains 33 clusters selected to evenly sample a range of 
X-ray luminosities and temperatures (2 to 9 keV), in- 
cluding a variety of dynamical states and core cooling 
ti mes. 

iCroston et al.1 (|2008| ) derived the radial gas density 
profiles for the REXCESS sample from X MM— Newton 
surface brightness profiles, using a non-par ametric depro- 
jection and PSF-deconvolution method. ICroston et al.l 

(2006) showed that this method accurately recovers the 
gas density profile observed in higher resolution Chandra 
data, down to 2'.'4 radius, for data of similar statisti- 
cal quality to REXCESS. Two clusters, RXCJ0956 and 
RXCJ2152, were excluded from the X-ray analysis be- 
cause they have multiple distinct components which pre- 
clude a one dimensional profile analysis, leaving a sample 
of 31 clusters. 

Relevant X-ray properties are listed in Table Q] The 
scaling radius rsoo (defined as the radius enclosing a 
mean overdensity of 500 times the critical density) and 
the corresponding cluster mass M c i was found from it- 
eratio ns about the M500 — Yx relation of Arnaud et al. 

(2007) . The cooling time was determined at a radius of 
0.03 r 50 o (about 10"-20"). The core gas density n e was 
measured at a radius of 0.008 r 50 o (about 3"-6") and 
was scaled by h(z)~ 2 to remove evolution effects. The 
dynamical state of the gas was characterized by (w), the 
standard deviation of the offsets w between the X-ray 
peak and the centroids of emission found for various radii 
( Bohringer et al. , i n prep). 

IPratt et al.1 (|2009f ) identified a subsample of REXCESS 
as cool core clusters using a cutoff in core gas density. 
We chose instead to use a cutoff in cooling time (at 
0.03 rsoo), since the cooling time can be compared to 
various dynamical time scales in the cluster. We defined 
cool core (CC) clusters to be those with t coo i < 2 Gyr (log 
tcooi < 9.3), marked with blue stars on the figures. Non- 
cool core (non-CC) c l usters have t coo i > 2 Gyr. Com- 
pared to IPratt et al.1 ((2009), this definition shifts only 
one cluster into the CC sample, namely RXCJ0211, see 

m 

HPratt et al.l (|2009l ) also identified clusters as morpho- 
logically disturbed if (w)> 0.01r 50 o. We use the same 
cut-off but with revised values of (if) from Bohringer et al. 
(in prep). Disturbed clusters are marked with red 
squares on the figures. In these clusters the gas distri- 
bution is less symmetric, likely due to a recent merger 
of sub-clusters and a younger dynamical age. In gen- 
eral, disturbed clusters do not have cool cores (although 
RXCJ1302 and RXCJ2319 are exceptions). 

For each cluster, the position of the X-ray peak emis- 
sion has been determined by Bohringer et al. (in prep), 
who describe the method in detail. Briefly, the peak po- 
sition for most clusters is the local maximum in a 4" 
smoothed image. In a few clusters the central surface 
brightness is flat, so the peak position is given as the 
center of a dipole fitted at 0.1 ^500. The positions were 
confirmed by visual inspection. These positions are more 
precise and acc urate than the low-resol ution ROSAT po- 
sitions listed in lBohringer et aTl (|2007l ). 

3. OPTICAL DATA 
3.1. Observations 



Table [5] lists our observations, which, with the excep- 
tion of two targets, were made with the 4.1 m South- 
ern Observatory for Astrophysical Research (SOAR) in 
Cerro Pachon, Chile. Most of the SOAR observations 
used the SOAR Opt ical Imager (SOL IWalker eTaIll2(M 
iSchwarz et al.l l2004[ ) . consisting of two thinned, back- 
illuminated 4096x4096 CCDs. The field of view is over 
5' across, with a binned pixel size of 0'.'154 from a natu- 
ral pixel size of 0'.'0767. The standard observing strategy 
was three 200 s exposures, dithered by 10" to span the 
gap between the two CCDs. The remainin g SOAR obser- 
vatio ns used the Goodman Spectrograph (i Clemens et al.l 
12001 in imaging mode, which utilizes a single Fairchild 
CCD. The field of view is 7'.2 in diameter, with a binned 
pixel size of 0'.'29 from a natural pixel size of 0'.'15. The 
standard observing strategy was a single 600 s exposure. 
One target (RXCJ1044) was observed at the 42 inch 
(1.1 m) John S. Hall Telescope at Lowell Observatory 
in New Mexico, using a 2048x2048 CCD camera with 
pixel size 0'.'59; we made 15 dithered exposures of 600 s 
each. All observations were made under clear photo- 
metric conditions; if the conditions or data were sus- 
pect in any way, the galaxy was reobserved. The seeing 
was typically 1". Finally, one target (RXCJ1311) was 
mistakenly omitted from observations but is one of the 
few to appear in the Sloan Dig ital Sky Survey (SDSS, 
lAdelman-McCarthv et "all [2008). SDSS data have suffi- 
cient resolution (pixel size 0'.'396) and depth for our mea- 
surements; a Luminous Red Galaxy spectral template 
was used to transform the r' data to the Cousins R filter 
(center ~630 nm) that we used for all other observations. 

The images were reduced in the standard fashion us- 
ing IRAF v.2.12.2, including bias subtraction, flat field 
correction, and cosmic ray removal. The dithered im- 
ages were aligned and median combined for each target, 
compensating for differing background levels if necessary. 
The background value of the combined image was mea- 
sured using APPHOT in a 1 to 2 arcminute radius an- 
nulus centered on the BCG with the radius chosen to 
avoid field objects. The data were corrected for airmass 
assuming an R band extinction coefficient of — 0.1, and 
calibr ated to the Vega magnitude scale u sing lLandoltl 
(1992) standard stars or lHamuv et al.l (|1992f ) spectropho- 
tometry standard stars. The data were converted to the 
AB scale with Rab = Rvega + 0.206; all magnitudes and 
colors are reported on the AB scale. 

The uncertainty in the final magnitudes is dominated 
by systematics. Because the BCGs are bright and cover 
many pixels, the random error on the summed flux and 
on the subtracted background is much smaller than the 
total systematic error in the flat-field corrections, stan- 
dard star calibration, airmass corrections, and variation 
in sky conditions throughout the night. We estimate 
the uncertainty due to these systematics to be around 
0.05 mag and assume this value for all observations. 

Images of a select subsample of the BCGs are 
shown in Fig ure [1] The entire sample is avail- 
able at http://www.calvin.edu/~dhaarsma/rexcess. 
The optical images are overlayed with contours of the 
XMM- Newton data, taken in the 0.5 - 2 keV band with 
all three detectors combined. The X-ray surface bright- 
ness images were corrected for vignetting and detector 
gaps, point sources were removed and smoothly refilled, 
and the images were smoothed with a 4" Gaussian. 
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3.2. Surface brightness profiles 

Figure shows the observed surface brightness profiles 
of the BCGs, corrected for Galactic extinction, to rest- 
frame, and for cosmological surface brightness dimming. 
T he slopes appear similar a cross the sample. 

iPostman &; Lauerl (|1995l ) characterized the profile 
using the slope of the integrated profile, a = 
d\og(L r ) / dlog(r) . They found that a depends on metric 
luminosity, and that this dependence could be removed 
to improve the utility of BCGs as standard candles. We 
measured a for this sample and found no correlation with 
metric luminosity, as well as an absence of low a val- 
ues. This confirms th e work of P. Lyman discussed in 
iBohringer et~aT1 (|200l . 

We also fit de Vaucouleurs functions to the profiles, but 
again found no correlation between profile shape (r e g) 
and metric luminosity. (A typical de Vaucouleurs func- 
tion is plotted in Figure [2]) Instead, we found the met- 
ric luminosity to be correlated with surface brightness 
at a given metric radius (the vertical offset in Figure [2]). 
Thus, these BCGs share the same profile shape but differ 
in their overall brightness. 

3.3. BCG identification 

Tables [3] and |4] list the optical properties of the BCGs. 
We classified the clusters into three categories. 1) Clus- 
ters where a central galaxy is clearly the brightest and 
has a generally smooth light distribution. 2) Clusters 
where a central galaxy is clearly the brightest, but has 
multiple light peaks appearing in the BCG envelope. 
These may be stars, foreground galaxies, infalling galax- 
ies on a radial path, member galaxies in low radius orbits, 
or mergers in progress. 3) Clusters in which there are two 
or more galaxies of similar brightness. Figure Q] includes 
examples of each type. 

We identified the BCG using aperture magnitudes, 
with two borderline cases decided by proximity to other 
galaxies or to the X-ray peak. The fields were inspected 
visually in the Digital Sky Survey out to 500 h" 1 kpc; if 
there was any ambiguity about the brightest galaxy or 
the X-ray peak position, the 5'-7' field of view was po- 
sitioned to cover other candidates. The brightness was 
measured for galaxies within this field of view (typically a 
200 h _1 kpc radius around the X-ray peak). We used an 
unmasked 50 h -1 kpc aperture in order to include light 
from the extended envelope and mergers in progress. One 
borderline case was RXCJ2234, in which three galax- 
ies had nearly the same brightness (within 0.01 mag) 
and none were near the X-ray peak, so we choose the 
one located in the middle of the group. The other 
case is RXCJ1311 (Abelll689), in which some galax- 
ies had a neighbor galaxy within the 50h _1 kpc aper- 
ture which boosted their brightness; measurements in a 
smaller (12 h -1 kpc) unmasked aperture identified two 
galaxies with similar brightness (within 0.06 mag) near 
the center of the group, so we choose the one located at 
the X-ray peak. 

BCGs located more than 0.03 rsoo from the X-ray peak 
are marked with an arrow in Figure [1] and the online 
images, and with a green X in Figures [3]|SJ Table 0] lists 
the 2MASS identification of the galaxy, but we were able 
to better resolve the crowded fields and determine a more 
accurate position of the BCG center, which we list in 



Table O 

3.4. Integrated Magnitudes 

Because of the crowded fields around BCGs, a simple 
aperture magnitude would be biased by extra flux from 
other galaxies in the field. Instead, we masked light from 
sources unrelated to the BCG, then fit elliptical isophotes 
to the remaining light distribution using the IRAF task 
ELLIPSE. We then integrated the isophotes numerically 
to get the total flux. 

Masked objects included stars, foreground galaxies, 
and cluster members not physically near the BCG. Some 
of the brightness peaks in the envelope, however, are the 
cores of cluster members in the process of merging with 
the BCG (such multiple nuclei are particularly visible 
in BCGs we classify as type 2). One approach would 
be to include the light from all objects that are already 
merging in the envelope. That is difficult, however, to 
implement, because it would require discriminating be- 
tween foreground galaxies and galaxies physically near 
the BCG, and determining a priori how much a galaxy 
must be apparently assimilated in order to include it with 
the BCG light. We chose instead to mask all brightness 
peaks except the BCG center, and to fit isophotes only 
to light in a smooth elliptical distribution around that 
center. This procedure provides a consistent measure 
without subjective decisions on which brightness peaks 
to include. As a test, we measured the magnitudes in 
a 50 h -1 kpc radius aperture without any masking, and 
found that the standard deviation of the absolute mag- 
nitudes increased from 0.29 to 0.34 mag; thus, the small 
masked aperture is a more consistent measure. 

When fitting ellipses, we fixed the central coordinates 
but allowed the ellipticity and position angle to vary with 
radius. If the ellipticity or position angle changed rapidly 
or the fit failed due to too much masking, we adjusted 
the masks to achieve a smooth fit. To prevent double 
counting we discarded any overlapping isophotes. 

The fitted isophotes were then integrated numerically, 
effectively replacing flux in the masked areas. During 
integration, we simulated a circular metric aperture by 
calculating the fraction of each ellipse inside the circular 
aperture. The apparent metric magnitude in a radius 
of 12/i _1 kpc (17 kpc for our assumed cosmology, about 
5"-15") is listed in Table [3] as R met . 

The apparent magnitude was converted to an absolute 
magnitude using the known redshi fts. We appl i ed the 
Galactic extinction corrections of Schlcgel ct al. (1998) 
as implemented in the calculator at the NASA Extra- 
galactic Database; corrections ranged from 0.04 to 0.40 
mag. We converted the magnitudes to rest-f rame using 
the kcorrect software (jBlanton fc Roweisl20 07?) , assuming 
a non-evolving Luminous Red Galaxy spectral template; 
corrections range from 0.05 to 0.18 mag. In Tables [3] 
and 2] we list magnitudes and colors with and without 
k-corrections so workers may apply their own SED and 
star formation history. 

3.5. Colors 

BCGs typically have similar colors because they have 
old stellar populations. To confirm this, we mea- 
sured R — K colors for the sample, making use of the 
Two Micron All Sky Survey Extended Source Catalog 
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(2MASSX. ISkrutskie et all [2006). In 2MASSX, the rec- 
ommended catalog magnitudes (Kk20fe) are isophotal 
(K = 20 mag arcsec" 2 with typical radii 10 to 50 kpc) in- 
side an elliptical aperture with bright neighbors masked 
out. To determine accurate colors, however, both bands 
must be measured in the same aperture with the same 
objects masked. We could cither use the R band mask 
and ellipses f ^3.4|) on the 2MASS image (giving a K to 
match Rmet), or use the 2MASSX mask and ellipses on 
the R band image (giving an Rk2q to match Kk20fe)- 
We used the latter method because of the difference in 
pixel scale between SOAR and 2MASS (0'.'154 vs. 1"). 
If the precise SOAR masks were applied to the low res- 
olution 2MASSX image, many more pixels of BCG light 
would be masked, leaving few pixels for the ellipse fit. In- 
stead, we masked the same objects in the SOAR image 
as were masked in 2MASSX and fixed the central coor- 
dinates, ellipticity, and position angle to the 2MASSX 
values during the ellip se fits. We then proceeded with 
integration as in £13.41 The resulting Rk20 values and 
R — K colors are listed in Table [4] 

Figures [3] and |4] show the R — K colors as they de- 
pend on absolute magnitude and on the J — K color 
from 2MASSX. With the exception of RXCJ1302 and 
RXCJ2319 (see , the colors have very little scatter 
and no trend with abolute magnitude. This consistency 
demonstrates that our methods have not introduced sig- 
nificant systematic errors relative to the 2MASS photom- 
etry. 

4. CORRELATION OF BCG AND GAS PROPERTIES 

4.1. Cluster mass vs. BCG mass 

Hierarchical formation models suggest that the mass 
of the BCG should be correl ated with th e mas s of 
its host cluster. For example, Wh ilev et al.1 (|2008f ) re- 
cently cal culated the expected depe ndence, using the 
models of iDe Lucia fc Blaizotl (|2007lh to be M B cg oc 
M° ; - 4 or M^ 5 depending on the feedback model used. 
A correlation between BCG luminosity and cluster X- 
ray l uminosity has long been noted in the literature 
(e.g. ISchombertiH9 88; Edge [l99lHEdge fc StewarflUMH: 
IHudson fe Ebelingj Il997ll Here we give an overview 
of recent results, which tend toward a shallower power 
law dependence than pre dicted by models. For exam- 
ple, I Whilev et al.l (|2008t ) measured the correlation as 



M BCG oc m; 



0.12±0.03 



for K band magnitudes inside a 



diameter of 37 kpc (radius of 13h 1 kpc). IBrough et al.l 
( 2008) saw a similar dependence of L BCG oc M c °- 11±01 ° 



at K band inside Uh" 1 kpc. lYang et al.l (|2008[) found 
Lbcg oc M°i 17 for galaxy groups. ILin fc Mohrl (|2004T l 
found a steeper dependence {Lbcg oc Af ° z ' 26±0 ' 04 ), but 
they used isophotal magnitudes rather than metric mag- 
nitudes. (Isophotal magnitudes favor larger, more mas- 
sive galaxies; in terms of Figure [21 isophotal magni- 
tudes make a horizontal cut through the profiles while 
metric magnitudes make a vertical cut.) Similarly, 
IPopesso et all (|2007D fo und L B cg p c M°/ 2 5 using r band 
Petrosion magnitudes. iStott et ail (|2008f) used K band 
isophotal magnitudes and found an even steeper rela- 
tion (— 1.1 ± 0.3 mag per decade of X-ray luminosity, or 
Lbcg oc £^ 44±0 ' 12 ) w hen fitting only to clusters with 
high X-ray luminosity. iMittal et al.l (|2009l ) used K band 



total magnitudes (extrapolating the Sersic profile) and 
found L B cg oc M°, 62±0 05 ; they note that total magni- 
tudes give a stronger correlation with cluster properties 
than smaller aperture magnitudes. 

For our sample, we measured Lbcg in R band inside 
YlhT 1 kpc. We confirm the weak trend between BCG 
luminosity and cluster mass. Figure [5] shows the trend, 
which has substantial scatter; Table [5] sho ws that the 
correl ation is marginally significant. Unlike IStott et"aL1 
(|2008h . we do not see the trend becoming more prominent 
at higher X-ray luminosities. Fitting in log-log space us- 
ing the bivariate correlated errors and intrinsic scatter 
(BCES(Y|X)) regression method (jAkritas fc Bershadvi 
11996ft . we find L B cg oc M cl 1%±om , consistent with the 
results found by others using metric magnitudes. 

4.2. Core gas density vs. BCG ellipticity 

We noticed a possible correlation between core gas den- 
sity and BCG ellipticity. Neglecting RXCJ2319 (see ©, 
the largest ellipticity among the cool core clusters is 0.29. 
The largest ellipticity among the non-CC clusters is 0.53, 
and half of the non-CC clusters have ellipticities larger 
than 0.29. This comparison suggests that BCGs in CC 
clusters are rounder than those in non-CC clusters. A K- 
S test shows that this trend is not statistically significant 
for the small numbers of our sample, so a larger sample 
is needed to test the trend suggested by these data. 

4.3. Core gas density vs. BCG location 

The offset between the X-ray peak and the BCG 
position is listed in Table [T] and plotted in Figure |H 
The XMM— Newton resolution prevents a precise mea- 
sure of this offset at the smallest separations, but 
XMM— Newton can clearly identify clusters where the 
BCG is located far from the X-ray peak. We find that 
three BCGs are located more than 0.035 rsoo (and more 
than 30") from the X-ray peak, all of which have multiple 
galaxies of similar brightness. That leaves 28/31 (90%) 
of the clusters with a galaxy closer than 0.035 rsoo- Note 
that the REXCESS sample is unbiased with respect to X- 
ray morphology, so this close agreement is not due to an 
overrepresentation of bright X-ray cores in the sample. 

Our high f r action of small offsets (90%) is higher than 
ILin fc Mohrl (|2004f ) found in a sample of 93 clusters: 
~65% within the same radius (0.02 -R200 in their Fig- 
ure 1). The difference from our results is likely because 
their sample is infrared-selected, has X-ray data drawn 
from several catalogs, and uses X-ray peak positions from 
lower resolution ROSA T data. 

We also found more small offsets than lLoubser et al.l 
(2009) found in a sample of 49 BCGs. They found 45% 
of BCGs falling within 20 kpc of the X-ray peak while we 
find 77% within that distance. This difference is likely 
due to different sample selection and to their exclusion 
of galaxies that fall near the peak if the peak is far from 
the literature X-ray position (their Table 7). 

Our sample i s mo re comparable to that of 
iSanderson etaTI (|200l . who study the LoCUSS 
sample of 65 X-ray selected clusters. They found 75% 
of BCGs falling within 0.04 rsoo of the X-ray centroid, 
while we found 90% within that distance of the X-ray 
peak. These percentages are consistent given Poisson 
statistics. 
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Our results agree with Hu dson et al.l (|2009f ). who study 
64 clusters in the HIFLUGCS sample. They found 88% 
to have a BCG within 50ft. -1 kpc of the X-ray peak, while 
we find 90% within that radius. 

Rec ent w ork bv IBildfell et all (12008ft . iRaffertv et all 
(|2008ft . and iSanderson et all (|2009h found that star for- 
mation in cool cores requires that the X-ray and galaxy 
centroids lie within ^20 kpc of each other. Our data 
agree; all 11 CC clusters have offsets less than 20 kpc, 
while 13 out of 20 non -CC have such small offsets. Like 
ISanderson et al.l ()2009f ) , we find that steep central X-ray 
profiles have small BCG offsets, similar to the correla- 
tion between central gas density with BCG offset shown 
i n Figure El 

IBildfell et all (|2008t ) find a correlation between BCG 
absolute B band magnitude and the offset from the X- 
ray peak, where BCGs are fainter when located far (up 
to 500 kpc) from X-ray peak. We see no such correlation 
at large offsets in R band. 

4.4. Gas dynamical state vs. BCG location 

We noticed a possible correlation between large values 
of (w) (indicating disturbed gas) and large offset between 
the BCG and the X-ray peak. RXCJ2048 and RXCJ2129 
are two of the three most disturbed clusters (the other is 
RXCJ2157), and are also two of the three largest BCG 
offsets (the other is RXCJ2234). The correlation is diffi- 
cult to quantify for small offsets, since many of our mea- 
sured offsets are consistent with zero. The correlation 
for large offsets is not surprising: clusters with disturbed 
gas have likely undergone recent mergers and the largest 
galaxies have not had time to settle into the core. 

4.5. Core gas density vs. BCG central stellar density in 
non-CC clusters 

4.5.1. Observed Trend 

We compared several optical properties of BCGs with 
various X-ray properties of the intracluster gas, and in 
most cases we saw little or no correlation. (Some of these 
are reported in Tableland describe above.) The tightest 
correlation, however, between X-ray and optical proper- 
ties was one we did not expect: between the core gas 
density of non-CC clusters and their aperture BCG lu- 
minosity, shown in Figure [7J 

The correlation appears to follow a power law, with 
outliers mainly toward higher core gas density. The 
rank order correlation statistics are reported in Ta- 
ble [5] for the non-CC sample (excluding RXCJ0211 and 
RXCJ2048, see Sj5j, the parameters show less than 0.16% 
chance of being uncorrelated, a detection stronger than 
3 sigma. A BCES (orthogonal) fit in log-log space finds 
n e cx £bcg° 4 - We also note that the correlation cannot 
be a spurious effect of the X-ray point spread function, 
since non-CC clusters have broad, flat centers that are 
well resolved by XMM— Newton. 

Note that the correlation is with aperture luminosity, 
not total luminosity. The aperture luminosity is pro- 
portional to the mass in the inner 12ft, -1 kpc, and so 
is a good proxy for the central stellar mass density of 
the galaxy (it is not a good proxy for the total mass 
of the galaxy). When discussing this trend, we will re- 
fer to Lbcg a s central stellar density. Thus, we discov- 
ered that, for non-CC clusters, the BCG central stellar 



density is larger in clusters with larger core gas densi- 
ties. We tested the correlation using several measures 
of BCG luminosity, but found that well-masked met- 
ric magnitudes give the tightest correlation (Table [5] 
shows P=0.0016). Well-masked isophotal magnitudes 
in R band also showed some correlation (P=0.015), 
but Kk20fe isophotal mags from 2MASSX showed lit- 
tle correlation (P=0.14). 2MASSX reports magnitudes 
with some bright objects masked, but not as completely 
masked as we could do with deeper, better-sampled im- 
ages in R band. 

We tested whether other related parameters would 
show the relationship more clearly, but found that core 
gas density and BCG metric luminosity had the strongest 
correlation (Table [5]). In place of BCG metric luminos- 
ity, we tried metric surface brightness. In place of core 
gas density, we tried cooling time and centroid shift. In 
no case was the cor relation as tight as the one shown. 
iCroston et al.l (|2008l ) reported the correlation of n e with 
X-ray centroid shift (w) (in opposing directions for the 
CC and non-CC samples), which in turn introduces a 
correlation of L bcg with (w) , but this is not as tight as 
Lbcg with n e . We tested if the correlation of n e and 
Lbcg could be tightened further by removing any addi- 
tional dependence on (w), but saw little improvement. 

We note that the correlation is driven by a pop- 
ulation of 5 clusters with the lowest core gas densi- 
ties: RXCJ0145, 0225, 2129, 2157, 2234. These 5 have 
not only the lowest core gas density, but are also the 
only clusters with multiple galaxies of similar bright- 
ness (which we classify as type 3). Three of the clusters 
(RXCJ0145, 2129, 2157) are among the most disturbed 
clusters in the sample. All three indicators suggest a 
young dynamic age, so this subsample points to cluster 
dynamics as key to the interpretation of the relation. 

4.5.2. Interpretation 

What could be the physical cause of this correlation? 
In no way do we suggest that the ICM gas density has 
a direct effect on the stars of the BCG; that would be 
unphysical. Rather, the explanation of the correlation 
must be based in some other physical process that affects 
both the core gas density and the central stellar density. 
What processes could be involved? 

BCGs and core gas are known to be related via the 
stars that form from the cooling gas. This causes BCGs 
in CC clusters to show more indicators of star for- 
mation than B CGs in non-CC cluste rs, including Ha 
emission (e.g. ICavagnolo et al.l [2008h . blue color (e.g. 
I Bildfell et al.ll2008ft. and younger stellar populations (e.g. 
lLoubser et al J 1 2 0*091) . Similar studies are being done for 
the REXCESS sample (Donahue et al. submitted) Our 
correlation, however, is seen in R band, which detects 
the mass of the old stellar population, not recent star 
formation. 

In fact, it seems unlikely that this old stellar popula- 
tion should be correlated with the presence or absence 
of a cool core, since the cool core gas and the stars 
have such different time scales. Prominent cool cores 
have radiative cooling times shorter than a billion years, 
whereas the vast majority of the stars formed several bil- 
lion years ago, based on studies of their color-magnitude 
relation and spectral energy distr i butio ns at z ~ 1 (e.g. 
lAndreon etahl 120081: IStott et al.l 12008ft . iCollins et all 
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(2009) suggest that even the assembly of BCGs is mostly 
complete before redshift 1. 

For these reasons, the correlation is not explained by 
the relationship between strong cool cores and star for- 
mation. Moreover, the correlation is in the non-cool core 
population, not the cool core population. In non-CC 
clusters, the gas is much too diffuse to cool and form 
stars. Some process other than star formation must be 
at work. 

The correlation is not due to overall cluster properties, 
such as total mass, X-ray luminosity, or gas temperature. 
It is not the case that more massive clusters simply have 
more luminous BCGs and denser core gas. While the 
BCG central stellar density does depend weakly on to- 
tal mass (Figure [5]) with a power of 0.18 ± 0.07, that 
trend cannot solely explain the much stronger depen- 
dence (power of 2.7 ±0.4) on core gas density. Moreover, 
the core gas density has no co rrelation with system tem- 
perature ((Croiton!eFld][200||) or with total cluster mass 
(Figure Table El . 

The correlation is not primarily due to the BCG grav- 
itational potential. It is true that the BCG is located 
near the X-ray peak in most cases (Figure El , and that 
a larger stellar density would increase the gravitational 
potential well. If this were the sole cause, however, the 
core gas density would be expected to be only modestly 
dependent on the central stellar density, but the observed 
relationship has a power of 2.7. In addition, two clusters 
with BCGs located far from the X-ray core also fall on the 
trend, implying that physical proximity is not required. 

The correlation might be due to a few clusters that 
have experienced a recent merger. As noted above, the 5 
clusters with lowest gas density include the three clusters 
with the most disturbed gas, suggesting a recent merger 
(Table [5] shows the strong correlation of gas density and 
(w)m nonCC clusters). These 5 clusters also have mul- 
tiple galaxies of similar brightness rather than a single 
dominant BCG, as expected for a recent cluster merger 
where the central galaxies have not yet merged into one 
BCG. After the merger, the resulting BCG would have 
higher total luminosity, but before the merger there is a 
time window when the gas density is low and the BCG 
total luminosity is low. This sounds like a promising ex- 
planation of the trend, but there is a problem. The Lbcg 
we report is a metric aperture luminosity in 12/i _1 kpc, 
not a total luminosity; as noted above, it is a better 
proxy for central stellar density than for total mass. How 
is central stellar density affected by a merger? Merg- 
ers between ellipticals are known to decrease the central 
density cusp, producing the well-known trend of lower 
central surface brightness in higher luminosity elliptical 
galaxies. Thus, the galaxies corresponding to low gas 
density are those with low stellar densities, which would 
be after their mergers rather than the cuspy galaxies seen 
in the window between gas merger and galaxy merger. 

We suggest that mergers are key to the explanation, 
but that it is the long-term merger history of the cluster 
rather than the most recent merger. Mergers of sub- 
clusters cause gas shocks which increase the gas entropy. 
More massive clusters have larger shocks and thus larger 
entropy gains. In the absence of galaxies and feedback, 
models show that the final entropy profile of the gas is 
determined by total cluster mass, leadi ng to a self- similar 
family of entropy profiles (reviewed in iVoitl [2005). The 



entropy at the core, however, tends to deviate from the 
self-similar relation. 

The core gas density can serve as a proxy for the devi- 
ation of the core entropy relative to other clusters of the 
same mass. In a self-similar cluster model, the ratio of 
core entropy Kq to the entropy K/\ at some large scale 
radius ta is 




The gas density n& at r& is the same for all clusters 
by definition. The ratio To/Ta can be affected by non- 
gravitational processes but is closely tied to gravitational 
potential, so To /Ta is not observed to vary much across 
self-similar clusters. Thus, in practice, the expression 
reduces to 



In other words, the central gas density is closely related 
to the relative enhancement of central entropy over what 
would be expected in a self-similar cluster, independent 
of cluster mass. A cluster whose central entropy was 
unusually large for its halo mass would therefore have an 
unusually low central density. 

Many physical processes can change the core entropy 
level, including radiative cooling, supernova or AGN 
feedback, and conduction. What we have discovered here 
is that, whatever the process, it also appears to set the 
central stellar density of the BCG. While most processes 
that affect gas entropy have little effect on stars, merg- 
ers affect both. The details of the merger history, such 
as the size of the subclusters, may be the cause of these 
deviations from the self-similar value. As stated above, 
mergers also affect stellar density, decreasing the cusp 
and reducing the central stellar density. Thus, merg- 
ers could increase gas entropy, decrease gas density, and 
decrease stellar density. In this way, the central stellar 
density would provide a record of the merger history of 
the cluster. The core gas in non-CC clusters has a very 
long cooling time and does not experience feedback, so it 
too preserves a memory of the shocks experienced in the 
merger history. In our sample, even disturbed clusters 
(red squares on the figures) fall on the trend, suggesting 
that the correlation is dominated by long-term merger 
history, not recent merger events. 

The CC sample does not follow a similar power-law 
trend. In these clusters, radiative cooling and feedback 
have a greater influence on the core gas, allowing the en- 
tropy to decrease and the gas density to increase beyond 
the level set by mergers. This increase does not have 
a particular correlation with BCG stellar density. The 
outliers on the upper side of the non-CC trend may be 
in transition to higher gas density via these processes. 

5. SPECIAL CASES 

RXC J0211: This cluster is classified by iPratt et al.l 
(2009) as non-CC based on its central gas density and 
central cooling time, but we classify it as CC based on 
its cooling time at 0.03 r$oo (see It has the lowest 
gas temperature of all the clusters in REXCESS, giving 
it a shorter cooling time than other clusters of similar gas 
density such as RXCJ0645. Like all CC clusters in REX- 
CESS, it hosts a central radio source, something that 
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some non-CC clusters lack (a full analysis of REXCESS 
radio properties will appear in Heidenreich et al, in prep). 
If RXCJ0211 is instead put in the non-CC sample, the 
trend between core gas density and BCG stellar density 
(Figure §4.5. 1|) is still present but a bit less significant 
(R=0.57, P=0.Q1). 

RXCJ0049 and RXCJ0225: In these clusters, the 
galaxy located at the X-ray peak is in the 2MASS Point 
Source catalog but not in the 2MASSX extended source 
catalog. Thus, we did not have the K band profile pa- 
rameters needed to measure the R — K color using the 
method in §3.51 

RXCJI302: This galaxy has a redder R — K color than 
the rest of the BCGs in the sample (Figure [3]) . Yet its 
Rmet abs magnitude (Figure [3]) and J — K color (Fig- 
ure [4j) are consistent with the rest of the sample. One 
possible explanation is dust internal t o the galaxy (found 
in other BCGs in CC clusters, e.g. lEgami et all 120061: 
lO'Dea et~a l. 2008). Assuming a typical extinction curve 
(|Gordon et alj 12003ft . about 0.5 mags of internal extinc- 
tion in R a bs,met would bring the R — K color in line with 
the rest of the sample while keeping J — K consistent 
with the sample. If such a change is made, the power 
fit in Figure [5] would change only slightly to 0.17±0.07. 
The fit in Figure [7J is unaffected because this cluster is 
not in the non-CC sample. 

RXCJ2048: This cluster is an outlier in the lower left of 
Figure [7] Of the clusters in the sample, it has the largest 
separation between the BCG position and X-ray peak 
(Figure H]). The BCG may be misidentified, although 
the only other candidate in the field is slightly fainter 
and located just as far from the X-ray peak on the oppo- 
site side of the X-ray centroid. Its X-ray properties are 
unique as well, having the most diffuse X-ray emission 
(Figure [T]) and the longest cooling time (Table [1} of the 
entire sample. It is possible the gas is not even in hydro- 
static equilibrium. Because of these multiple issues, we 
leave it out of the fit in Figure [7J and out of the statistics 
in Table [5j 

RXCJ23I9: This galaxy has a much fainter absolute 
magnitude than the rest of the BCGs in the sample (Ta- 
ble [3]) . We reobserved the galaxy on another night and 
confirmed the faint R magnitude. A misidentification 
is unlikely, since the galaxy is clearly the brightest in 
the field and coincident with the X-ray peak. A faint 
magnitude and red color could be caused by an incor- 
rect redshift , but the optical redshift z = 0.0984 from 
IGuzzo et al.l ([2009) is confirmed by our own observations 
of faint emission lines in the BCG at z = 0.0979 (Don- 
ahue et al. submitted). In addition, its R — K and 
J — K colors are extremely red (Figures [3l [4|, while its 
J — H color from 2MASSX is consistent with the sample. 
That suggests an excess of K emission, per haps similar 
to th e obscured Seyfert BCG in Abell 1068 ([Edge et al.l 
2002), except the optical spectrum shows very typical 
line widths rather than the strong lines of AGN activity. 
Thus, we do not have a good explanation for this object, 
but because of its puzzling nature we leave it out of the 
fits in Figures [3] and [7J and out of the statistics in Table[SJ 

6. SUMMARY 

Using the REXCESS sample of galaxy clusters, we 
investigated the relationships between the BCG stellar 



mass and the properties of the X-ray emitting gas. We 
confirmed the weak correlation seen by others between 
the BCG luminosity and the total cluster mass, and the 
close proximity of the BCG to the X-ray peak. 

We detected a trend among the non-CC clusters in 
which the core gas density increases with the BCG aper- 
ture luminosity, a proxy for the BCG central stellar mass 
density. This trend is much clearer when the BCG lu- 
minosity is determined from well-masked aperture mag- 
nitudes rather than from unmasked or isophotal mag- 
nitudes. This trend holds even in cases where the gas 
is disturbed or the BCG is located far from the central 
region. We argue that the core gas density is an indi- 
cator of the deviation of the central entropy from the 
self-similar value. Thus, the core gas entropy and the 
central BCG stellar density appear to be more closely 
related than previously thought. We suggest that clus- 
ter mergers could be the underlying cause, since mergers 
can both increase gas entropy and decrease BCG central 
stellar density. If so, this trend could set important con- 
straints on models of the central gas in clusters and of 
BCG formation. 

These results are based on a sample of only 31 clus- 
ters, and need to be co nfirmed using other larg er cluster 
samples. For example, iCavagnolo et al.l (|2009f) recently 
fit entropy profiles to 239 clusters from the Chandra X- 
ray Observatory archive. The optical data corresponding 
to this or other large samples of X-ray selected clusters 
would allow an investigation of the trend with a statisti- 
cally significant sample size. If the trend is real and our 
interpretation is correct, it should be present indepen- 
dent of how the clusters are selected. 
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TABLE 1 
X-ray Properties 



Cluster 
(1) 



Alt name 



(2) 



(3) 



^SOO 
(kpc) 
(4) 



n e h(z) 2 
(cm -3 ) 
(5) 



log(i coo ;) log({«;)) log(M c j) X peak coordinates Offset 

(log(yr)) (log(rsoo)) (log(M Q )) J2000 (kpc) 

(6) (7) (8) (9) (10) 



RXCJ0003. 
RXCJ0006. 
RXCJ0020. 
RXCJ0049. 
RXCJ0145. 
RXCJ0211. 
RXCJ0225. 
RXCJ0345. 
RXCJ0547. 
RXCJ0605. 
RXCJ0616. 
RXCJ0645. 
RXCJ0821. 
RXCJ0958. 
RXCJ1044. 
RXCJ1141. 
RXCJ1236. 
RXCJ1302. 
RXCJ1311. 
RXCJ1516. 
RXCJ1516. 
RXCJ2014. 
RXCJ2023. 
RXCJ2129. 
RXCJ2149. 
RXCJ2157. 
RXCJ2217. 
RXCJ2218. 
RXCJ2234. 
RXCJ2319. 



8+0203 
0-3443 
7-2542 
4-2931 

0- 5300 
4-4017 

1- 2928 

7- 4112 

6- 3152 

8- 3518 
8-4748 

4- 5413 
8+0112 

3- 1103 

5- 0704 

4- 1216 

7- 3354 

8- 0230 

4- 0120 
3+0005 

5- 0056 
8-2430 

0- 2056 
8-5048 

1- 3041 

4- 0747 
7-3543 

6- 3853 

5- 3744 

6- 7313 



A2700 
A2721 
A0022 
S0084 
A2941 
A2948 

S0384 
A3364 
A3378 

A3404 
A0653 
A0907 
A1084 
A1348 
A0700 
A1663 
A1689 
A2050 
A2051 

S0868 
A3771 
A3814 
A2399 
A3854 
A3856 
A3888 
A3992 



0.0924 
0.1147 
0.1410 
0.1084 
0.1168 
0.1008 
0.0604 
0.0603 
0.1483 
0.1392 
0.1164 
0.1644 
0.0822 
0.1669 
0.1342 
0.1195 
0.0796 
0.0847 
0.1832 
0.1181 
0.1198 
0.1538 
0.0564 
0.0796 
0.1184 
0.0579 
0.1486 
0.1411 
0.1510 
0.0984 



8761! 
105911° 

10451! 
807±7 

10891! 
685±| 
693±g 
688±4 

11331! 

1045l| 
939l| 

12791? 
755l| 
10771" 
931I2 
885I2 
7531! 
8421^ 

131911 
989l| 
9271! 

115511 
7391! 
900lg 
88611 
75111 

102211 

H30lg 

12831s 
78815 



0.0162+0. 
0.0079+0. 
0.0085+0. 
0.0182+0. 
0.0038+0. 
0.0214+0. 
0.0038+0. 
0.0547+0. 
0.0088+0. 
0.0739+0. 
0.0119+0. 
0.0245+0. 
0.0133+0. 
0.0503+0. 
0.1015+0. 
0.0701+0. 
0.0125+0. 
0.0347+0. 
0.0465+0. 
0.0109+0. 
0.0104+0. 
0.1291+0. 
0.0092+0. 
0.0052+0. 
0.0549+0. 
0.0034+0. 
0.0182+0. 
0.0126+0. 
0.0063+0. 
0.0571+0. 



0007 
0008 
0007 
0014 
0005 
0008 
0005 
0010 
0005 
0013 
0006 
0014 
0008 
0023 
0018 
0010 
0007 
0006 
0011 
0006 
0008 
0023 
0009 
0005 
0011 
0003 
0009 
0010 
0005 
0018 



9.47+0.02 
9.82+0.02 
9.67+0.02 
9.37+0.02 
9.92+0.02 
9.23+0.02 
9.60+0.02 
9.07+0.02 
9.65+0.02 
8.95+0.02 
9.64+0.02 
9.37+0.02 
9.60+0.02 
9.08+0.02 
8.87+0.02 
8.94+0.02 
9.41+0.02 
9.20+0.02 
9.16+0.02 
9.61+0.02 
9.66+0.02 
8.74+0.02 
9.61+0.02 
9.91+0.02 
8.92+0.02 
10.02+0.02 
9.45+0.02 
9.50+0.02 
9.80+0.02 
8.80+0.02 



-2.49+0.12 
-1.89+0.05 
-2.20+0.05 
-2.64+0.13 
-1.52+0.02 
-2.34+0.09 
-1.92+0.05 
-2.28+0.07 
-2.15+0.03 
-2.23+0.03 
-1.88+0.05 
-2.41+0.04 
-2.35+0.12 
-2.47+0.08 
-2.14+0.02 
-2.27+0.05 
-2.28+0.05 
-1.82+0.02 
-2.40+0.03 
-2.43+0.05 
-1.75+0.03 
-2.24+0.02 
-1.78+0.04 
-1.38+0.17 
-2.47+0.06 
-0.97+0.97 
-2.74+0.49 
-1.81+0.01 
-2.12+0.03 
-1.66+0.02 



14.321 



+0.008 



14.5773 
14.57lj 



-n.oo* 

7+0.012 
-0.013 
,+0.007 
-0.007 

14 221+ ' 012 
L1.ZZL_ Q q 12 

14 614+ ' 008 

14 003+ - 007 
,+0.014 



0.015 



14.003" 
13 992+ 011 



14.680i°;°° 7 7 

1+0.007 
-0.007 
14.421+ - 008 



14.57r 



-0.00K 
,+0.008 
-0.008 
,+0.010 
-0.011 
,+0.022 
-0.020 
5+0.004 
-0.004 
-+0.004 
-0.004 
i+0.011 
-0.002 
1-0.007 
-0.007 

14.8931!;!!! 

-,+0.005 
-0.005 
-+0.009 
-0.008 
1-0.005 
-0.005 



14.846Z 
14.123 j 
14.622^ 
14.419j 
14.345:! 
14.118 j 
14.265"* 



14.490 
14.405 
14.707 
14.084l!;!H 
14.351 



14.347 



+0.011 
0.012 

+0.006 



14.106"* 



-0.00<> 
,+0.009 
3 -0.009 
14 ^4fi+°-006 

14 Q4Q+0.005 
14.843_ a005 

,+0.009 



14.186" 



0.009 



00:03:49.7 
00:05:59.9 
00:20:42.2 
00:49:23.0 
01:44:59.7 
02:11:24.8 
02:25:09.3 
03:45:46.2 
05:47:38.4 
06:05:54.2 
06:16:51.7 
06:45:29.3 
08:21:50.9 
09:58:22.3 
10:44:33.0 
11:41:24.4 
12:36:41.3 
13:02:53.3 
13:11:29.5 
15:16:18.1 
15:16:44.2 
20:14:51.7 
20:22:58.8 
21:29:40.9 
21:49:07.6 
21:57:29.5 
22:17:45.5 
22:18:40.3 
22:34:27.1 
23:19:40.2 



+02:03:58 
-34:43:23 
-25:42:25 
-29:31:14 
-53:01:03 
-40:17:28 
-29:28:36 
-41:12:14 
-31:52:12 
-35:18:09 
-47:47:40 
-54:13:40 
+01:11:52 
-11:03:54 
-07:04:09 
-12:16:37 
-33:55:37 
-02:31:00 
-01:20:28 
+00:05:28 
-00:58:12 
-24:30:20 
-20:56:56 
-50:48:55 
-30:42:05 
-07:47:55 
-35:43:30 
-38:54:06 
-37:44:02 
-73:13:38 



5+ 6 
14+ 8 
32+ 9 
5+ 7 
23+ 8 
1+ 7 
4+ 4 
4+ 4 
39+10 
7+ 9 
14+ 8 
10+11 
6+ 6 
17+11 
4+ 9 
5+ 8 
8+ 6 
17+ 6 
1+12 
16+ 8 
6+ 8 
5+10 
4+ 4 
51+ 6 
4+ 8 
12+ 4 
10+10 
29+ 9 
112+10 
5+ 7 



Note. - (1) Name of cluster in REXCESS catalog. (2) Alternate name of cluster. (3) Redshift of the cluster 



Croston et al 



Croston et al 



(4) Th e radius of the cluster enclosing a mean overden s ity of 500 times the critical density, found by 
(]2008j ) using the M 500 - Y x relation of lArnaud et al.l (|2007l) . (5) Gas density at 0.008 r 500 , from 
(2008). (6) Gas cooling time at 0.03 r 500 , from iCroston et all (|2008l ): we classify cool cores as t coo i < 2 x 10 9 yr or 
log(t coo ;) < 9.3. (7) Standard deviation of centroid s hifts, from Bohring er et al. (in prep); we classify dist urbed clusters 
as (w) > 0.01 r5oo- (8) Total cluster mass, found bv lPratt et all ([2009) using the M500 — Yx relation of lArnaud et al.l 
(120071) . (9) Coordinates of the peak of the X -ray emission, from Bohringer et al. in prep, with an uncertainty of 4". 
(10) Offset of the selected BCG from the peak X-ray position. 
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TABLE 2 

Observations 



Date 



Targets 



Telescope Instrument 



2007 September 15 

2007 October 11 

2008 March 8 
2008 July 6 
2008 July 7 
2008 October 4 

2008 November 2 

2009 January 15-16 
2009 April 17 
2009 April 27 



RXCJ 2217 
RXCJ 0345 
RXCJ 0821 
RXCJ 2014 
RXCJ 0006 
RXCJ 2149 
RXCJ 0003 
RXCJ 1044 
RXCJ 0616, 1141 
RXCJ 1236, 1302 
RXCJ 1311 



2218, 2234, 2319 
0547, 0605, 0645 
0958, 1516.3, 1516.5 
2023, 2048, 2129, 2157 
0020, 0049, 0145, 0211, 0225 



SOAR 
SOAR 
SOAR 
SOAR 
SOAR 
SOAR 
SOAR 
Lowell 
SOAR 
SOAR 
SDSS 



SOI 
SOI 
SOI 
SOI 
SOI 

Goodman 
Goodman 
Hall 42" 
Goodman 
Goodman 



TABLE 3 
Optical Properties ■ 



metric: 



BCG coord R met R me t,abs L B CG PA 

J2000 ( ma g) observed rcstframe (log(Z,Q)) (deg) 

(1) (2) (3) (4) (5) (6) 



(7) 



(mag arcsec 
(8) 



00:03:49.6 


+02:04:00 


15.06 


-23.07 


-23.16 


11.06+0.02 


-81 


0.45 


21 


.24 


00:05:59.6 


-34:43:17 


15.59 


-23.04 


-23.16 


11.06+0.02 


-77 


0.51 


21 


.18 


00:20:43.1 


-25:42:28 


16.23 


-22.89 


-23.04 


11.01+0.02 


24 


0.47 


21 


22 


00:49:22.8 


-29:31:12 


15.53 


-22.98 


-23.09 


11.03+0.02 


-25 


0.04 


21 


.92 


01:44:58.9 


-53:01:12 


16.22 


-22.46 


-22.58 


10.83+0.02 


74 


0.21 


22 


.15 


02:11:24.8 


-40:17:28 


15.64 


-22.70 


-22.80 


10.91+0.02 


-62 


0.23 


22 


.09 


02:25:09.0 


-29:28:38 


14.61 


-22.55 


-22.61 


10.84+0.02 


-76 


0.31 


22 


.16 


03:45:46.0 


-41:12:16 


13.87 


-23.28 


-23.34 


11.13+0.02 


52 


0.09 


21 


.77 


05:47:37.7 


-31:52:24 


16.36 


-22.88 


-23.05 


11.01+0.02 


19 


0.09 


22 


.05 


06:05:53.9 


-35:18:08 


16.08 


-23.01 


-23.16 


11.06+0.02 


-73 


0.11 


21 


.93 


06:16:51.7 


-47:47:45 


15.22 


-23.45 


-23.57 


11.22+0.02 


83 


0.53 


20 


.73 


06:45:29.5 


-54:13:37 


16.26 


-23.23 


-23.41 


11.16+0.02 


65 


0.31 


21 


.19 


08:21:50.7 


+01:11:49 


15.21 


-22.66 


-22.74 


10.89+0.02 


-18 


0.34 


22 


.30 


09:58:22.0 


-11:03:51 


16.56 


-22.96 


-23.14 


11.05+0.02 


-27 


0.26 


21 


.60 


10:44:32.9 


-07:04:07 


16.39 


-22.62 


-22.76 


10.90+0.02 


2 


0.29 


22 


.07 


11:41:24.2 


-12:16:37 


15.49 


-23.24 


-23.37 


11.14+0.02 


-7 


0.27 


21 


.40 


12:36:41.3 


-33:55:32 


14.91 


-22.87 


-22.95 


10.97+0.02 


-16 


0.13 


22 


.14 


13:02:52.6 


-02:30:59 


15.15 


-22.78 


-22.86 


10.94+0.02 


78 


0.22 


22 


16 


13:11:29.5 


-01:20:28 


16.63 


-23.11 


-23.31 


11.12+0.02 


58 


0.13 


21 


.49 


15:16:17.9 


+00:05:21 


15.87 


-22.83 


-22.95 


10.98+0.02 


43 


0.41 


21 


.57 


15:16:44.2 


-00:58:09 


15.80 


-22.94 


-23.06 


11.02+0.02 


-45 


0.29 


21 


.59 


20:14:51.7 


-24:30:22 


16.10 


-23.33 


-23.51 


11.20+0.02 


37 


0.21 


21 


.38 


20:22:59.1 


-20:56:56 


14.14 


-22.86 


-22.91 


10.96+0.02 


-66 


0.09 


21 


.95 


20:48:11.6 


-17:49:03 


15.83 


-23.40 


-23.56 


11.22+0.02 


16 


0.14 


21 


.45 


21:29:42.4 


-50:49:26 


15.10 


-22.69 


-22.76 


10.90+0.02 


10 


0.19 


22 


.20 


21:49:07.4 


-30:42:05 


15.29 


-23.42 


-23.55 


11.21+0.02 





0.19 


21 


.36 


21:57:29.4 


-07:47:44 


14.45 


-22.61 


-22.67 


10.86+0.02 


80 


0.42 


22 


.55 


22:17:45.8 


-35:43:29 


16.00 


-23.25 


-23.41 


11.16+0.02 


-22 


0.20 


21 


.11 


22:18:39.4 


-38:54:02 


16.03 


-23.09 


-23.24 


11.09+0.02 


-35 


0.36 


21 


.24 


22:34:24.6 


-37:43:31 


16.47 


-22.81 


-22.98 


10.98+0.02 


18 


0.15 


21 


.84 


23:19:40.5 


-73:13:36 


16.39 


-21.88 


-21.99 


10.59+0.02 





0.26 


22. 


.72 



Note. — All magnitudes and colors are on the AB scale and corrected for Galactic dust extinction. (1) 
BCG coordinates (2) R band apparent magnitude in metric radius of 12 h — 1 kpc, uncertainty of 0.05 mag. 
(3,4) R band absolute magnitude in metric radius of 12 h — 1 kpc, uncertainty of 0.05 mag. Column 4 is 
converted to restframc. (5) Luminosity of BCG corresponding to Rmet,abs converted to rcstframe. (6,7,8) 
Position angle (degrees E of N), ellipticity, and corrected surface brightness of the isophote with semi-major 
axis 12 h~ 1 kpc. 
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TABLE 4 
Optical Properties - Isophotal 



2MASS Name 


Class 


Apcr 


Rk20 


R — 


K color 


J — K color 


SMA 










(") 


(mag) 


observed 


restframe 


restframe 


(kpc) 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2MASX J00034964+0203594 


1 


20.2 


14.57 


1.31 


1.08±0.10 


0.18±0.10 


40 


2MASX J00055975- 


-3443171 


1 


12.4 


15.46 


1.18 


0.90±0.10 


-0.05±0.10 


41 


2MASX J00204314- 


-2542284 


2 


18.7 


15.32 


1.39 


1.01±0.09 


0.14±0.09 


47 


2MASS J00492286- 


-2931124 


2 












30 


2MASX J01445891- 


-5301110 


3 


10.6 


16.01 


1.17 


0.88±0.13 


0.15±0.16 


25 


2MASX J02112484- 


-4017261 


1 


15.3 


15.33 


1.23 


0.98±0.11 


-0.03±0.12 


28 


2MASS J02250904- 


-2928383 


3 












29 


2MASX 


J03454640- 


-4112149 


1 


22.3 


13.61 


1.18 


1.03±0.07 


0.08±0.06 


28 


2MASX 


J05473773- 


-3152237 


1 


11.5 


16.00 


1.34 


0.94±0.16 


0.11±0.18 


26 


2MASX 


J06055401- 


-3518081 


1 


9.8 


15.87 


1.44 


1.07±0.09 


0.21±0.11 


25 


2MASX 


J06165166- 


-4747434 


1 


18.0 


14.74 


1.16 


0.87±0.10 


0.04±0.11 


65 


2MASX 


J06452948- 


-5413365 


1 


14.7 


15.73 


1.33 


0.87±0.13 


0.13±0.15 


39 


2MASX 


J08215065+0111495 


1 


12.9 


15.06 


1.22 


1.02±0.11 


0.10±0.12 


21 


2MASX 


J09582201- 


-1103500 


1 


11.3 


16.03 


1.56 


1.10±0.12 


0.13±0.13 


10 


2MASX 


J10443287- 


-0704074 


1 


7.1 


16.49 


1.26 


0.90±0.12 


0.04±0.13 


25 


2MASX 


J11412420- 


-1216386 


1 


11.2 


15.37 


1.34 


1.03±0.09 


0.16±0.09 


38 


2MASX 


J12364125- 


-3355321 


1 


16.1 


14.62 


0.90 


0.71±0.09 


-0.06±0.09 


23 


2MASX 


J13025254- 


-0230590 


1 


14.4 


14.97 


1.61 


1.40±0.09 


0.13±0.08 


24 


2MASX 


J13112952- 


-0120280 


1 


8.2 


16.32 


1.48 


0.96±0.12 


0.13±0.13 


26 


2MASX 


J15161794+0005203 


1 


10.1 


15.76 


1.36 


1.05±0.11 


0.14±0.12 


34 


2MASX 


J15164416- 


-0058096 


1 


11.6 


15.60 


1.17 


0.86±0.11 


0.02±0.12 


33 


2MASX 


J20145171- 


-2430229 


1 


9.4 


15.84 


1.28 


0.86±0.10 


0.06±0.11 


36 


2MASX 


J20225911- 


-2056561 


2 


19.9 


13.77 


0.79 


0.65±0.09 


-0.09±0.08 


25 


2MASX 


J20481162- 


-1749034 


3 


9.3 


15.63 


1.27 


0.87±0.10 


O.OliO.ll 


36 


2MASX 


J21294244- 


-5049260 


3 


11.6 


15.27 


1.03 


0.84±0.10 


-0.08±0.11 


25 


2MASX 


J21490737- 


-3042043 


1 


14.0 


14.91 


1.12 


0.81±0.09 


0.04±0.08 


38 


2MASX 


J21572939- 


-0747443 


3 


10.3 


14.72 


1.13 


0.99±0.07 


0.01±0.06 


18 


2MASX 


J22174585- 


-3543293 


1 


13.6 


15.25 


1.32 


0.91±0.10 


0.09±0.11 


15 


2MASX 


J22183938- 


-3854018 


2 


14.1 


15.50 


1.23 


0.85±0.09 


0.08±0.09 


51 


2MASX 


J22342463- 


-3743304 


3 


11.0 


15.88 


1.24 


0.83±0.10 


O.OOiO.ll 


33 


2MASX 


J23194046- 


-7313366 


1 


9.7 


16.44 


1.54 


1.30±0.13 


0.43±0.17 


16 



Note. — All magnitudes and colors (including 2MASS) arc on the AB scale and corrected for Galactic 
dust extinction. (1) 2MASS name of selected BCG (sec J5] for RXCJ0049 and RXCJ0225). (2) BCG 
classification (sec i|3.3l (3) Semi-major axis of aperture used for Ft — K color, i.e. the 2MASSX isophotc 
where K — 20 magarcscc -2 . (4) R band apparent magnitude in the 2MASSX K20 aperture, with masking to 
match 2M ASSX, uncertainty 0.05 mag. (5,6) R — K color, both bands measured in same aperture. Column 
6 is converted to restframe. (7) 2MASS J — K color, converted to restframe. (8) Semi-major axis of the 
isophote at 22.5 magarcsec~ 2 . 



TABLE 5 

Correlation statistics 



Parameters 


CC Clusters 


Non-CC Clusters 


All Clusters 




R 


P 


R 


P 


R 


P 


Optical - Optical 














if, - K vs. LsCG 


0.41 


0.24 


0.12 


0.64 


0.13 


0.49 


R - K vs. J - K 


0.73 


0.016 


0.70 


0.00078 


0.70 


0.000020 


X-ray - X-ray 














n e vs. M c i 


0.37 


0.29 


0.095 


0.70 


0.067 


0.73 


n e vs. (to) 


0.33 


0.35 


-0.74 


0.00026 


-0.56 


0.0017 


X-ray - Optical 














M cl vs. Lbcg 


0.16 


0.65 


0.46 


0.047 


0.38 


0.040 


n e vs. e 


-0.091 


0.80 


-0.067 


0.78 


-0.26 


0.17 


(w)vs. Lbcg 


-0.37 


0.29 


-0.44 


0.060 


-0.46 


0.012 


n e vs. Lbcg 


0.36 


0.31 


0.67 


0.0016 


0.56 


0.0017 



Note. — R is the Spearman rank-order correlation coefficient. P is the probability 
that the parameters are not correlated. Correlations with the offset between the BCG 
and X-ray peak are not listed because many offsets arc consistent with zero. Clusters 
RXCJ2319 and RXCJ2048 have been left out of all statistics, sec J5] 
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Fig. 1. — Images of REXCESS BCGs. The complete sample can be viewed at http://www.calvin.edu/~dhaarsma/rexcess/ The four 
clusters shown here illustrate our classification scheme: RXCJ0211 is type 1 (brightest, smooth), RXCJJ0049 is type 2 (brightest, lumpy), 
RXCJ2157 is type 3 (multiple galaxies of similar brightness) with the X-ray peak clearly selecting one galaxy, and RXCJ2048 is type 3 with 
the X-ray peak far from any galaxy. North is up and east left. Greyscale: R band. Vertical bar: 50 kpc at redshift of cluster. Contours: 
X-ray surface brightness from XMM— Newton, 0.5 — 2 keV, with contours increasing in steps of y2. 
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Fig. 2. — R band surface brightness profiles of the BCGs, including k-corrections, extinction corrections, and corrections for surface 
brightness dimming. The data are plotted only for radii greater than 3 arcsec and less than twice the isophotal radius (isophote of 22.5 
mag arcsec -2 ). Solid lines are the cool core population, dotted lines are the remaining BCGs. The vertical dashed line indicates the metric 
radius of 12/i~ 1 kpc. The thick line in the upper right indicates the shape of a typical de Vaucouleurs profile (effective radius of 60 kpc). In 
the online version of the paper, the colors of the lines indicate the metric luminosity of the BCG, showing the correlation with the surface 
brightness at the metric radius. 



1 .6 
1 .4 

1 .2 
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Fig. 3. — BCG color vs. absolute metric magnitude in AB units, restframe, corrected for Galactic dust extinction. Symbols indicate 
special sub-samples: cooling time less than 2 Gyr (blue star), disturbed X-ray emission (red square), and large separation between BCG 
and X-ray peak (green X). In general, the BCGs have similar color and no dependence of color on magnitude; the exceptions are RXCJ1302 
and RXCJ2319, see (0 
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Fig. 4. — BCG colors in AB units, restframe, corrected for Galactic dust extinction. Symbols same as Figure[3] The dotted lines indicate 
reddening due to dust internal to the BCG; the length of the line indicates extinction from Ay = to 1 for the outliers RXCJ1302 and 
RXCJ2319, see ij5] 




Fig. 5. — The BCG metric luminosity (radius 12/i 1 kpc) depends slightly on total cluster mass. Symbols same as Figure [3] The line 
shows the BCES(Y|X) fit to the data (excluding RXCJ2319, see SJ5) with a power law of 0.18 ± 0.07. 
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Fig. 6. — In cool core clusters, the BCG is co-located with the peak of the X-ray emission; 90% of BCGs in the sample are located within 
0.035 rsoo- Symbols same as in Figure [3] 
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Fig. 7. — The gas density of the cluster at 0.008 rsoo vs. the BCG metric luminosity (effectively the stellar mass density inside a radius 
of 12/i -1 kpc). Symbols same as Figure[3] The line is a BCES(orthogonal) fit to the non-cool core population (excluding RXCJ2048, sec 
tJSJ , with a power law of 2.7 ± 0.4. 



16 



Haarsma et al. 



0.100 



OT 0.010 
o 



ft • 



♦ 



SI 



0.001 







14 



,14 



3x10' 

Total Cluster Mass (M ) 







15 



Fig. 8. — The gas density of the cluster at 0.008 rsoo vs. the total cluster mass. Symbols same as Figure [3] There is no correlation 
between the core gas density and the overall cluster mass. 



